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We use exact diagonalization to determine the spectrum of reduced Hamiltonians based on renor- 
malization group flows to strong coupling. For the half-fiUed two-leg Hubbard ladder we reproduce 
the known insulating d-Mott groundstate with spin and charge gaps. For the saddle point regions of 
the two-dimensional Hubbard model near half-filling we find a crossover to a similar strong coupling 
state, which truncates the Fermi surface near the saddle points. At lower scales d-wave supercon- 
ductivity appears on the remaining Fermi surface. 
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In recent years the application of renormalization 
group (RG) methods to interacting fermions in low di- 
mensions has made great progress. In one dimension 
(ID) where the Fermi surface (FS) is a discrete set of 
points, the method has a long history going back to the 
pioneering work of Menyhard and Solyom jj]. In two 
dimensions (2D) the FS is a continuous curve and the 
discrete set of coupling constants is replaced by a 4- 
point vertex which is a function of three independent 
momenta. Zanchi and Schulz ^] developed a numerical 
method to evaluate the functional renormalization group 
flow of the scattering vertex based on a break-up of the 
FS into N patches. The resulting coupled differential 
equations can be integrated numerically to obtain the 
RG flow. This method has been applied to the Hub- 
bard model on a square lattice with nearest-neighbor 
(n.n.) and next-nearest neighbor (n.n.n.) hopping us- 
ing a Wilsonian scheme which successively integrates out 
high energy modes down to an infrared cutoff |^ |^ Q . 
The RG method has the advantage that it allows an unbi- 
ased treatment of competing instabilities including even 
Stoner ferromagnetism, when the recent reformulation 5], 
with the temperature as the flow parameter, is used. 

A serious limitation to RG schemes occurs when the 
scattering vertices flow to strong coupling as the RG scale 
is lowered. Although the one-loop flow may still be a 
good approximation even when the interaction strength 
has become comparable to the kinetic energy, self energy 
corrections A] typically become important shortly after 
that and higher order terms or non-perturbative effects 
might not be negligible any longer. Thus these RG flows 
cannot be extrapolated to scales below the divergence. 
Generally those are the most interesting cases as diver- 
gent RG flows signal a breakdown of the perturbative 
normal state. In the special case of ID, bosonization al- 
lows a complete analysis of the strong coupling regime. 
In higher dimensions, in cases where a single instabil- 
ity dominates the RG flow, a mean fleld approximation 
can be introduced to describe the strong coupling phase. 
This can be corrected for fluctuations around the mean 



fleld to obtain a satisfactory description of the ordered 
phase. Yet in some especially interesting cases, there can 
be several mutually reinforcing instabilities so that it is 
not clear how to formulate a mean field approximation. 

In this Letter we propose a new and more general 
approach to determine the strong coupling phase that 
emerges on scales below divergent RG flows. The start- 
ing point is to recognize that a mean fleld approxima- 
tion based on factorization of the interaction term in the 
Hamiltonian, can also be expressed as an exact solution of 
a reduced Hamiltonian e.g. BCS theory was formulated 
as an exact solution for a reduced Hamiltonian in which 
only terms describing the scattering of Cooper pairs with 
zero total momentum were kept 6]. Note these are the 
most divergent terms in the RG flow of the 4-point ver- 
tex. Thus we can regard the mean field approximation 
as contained in the solution of a reduced Hamiltonian 
whose form is derived from the diverging RG flow. In 
those cases mentioned above the reduced Hamiltonian is 
more complex and cannot be simply solved analytically. 
However it is amenable to numerical diagonalization and 
this is the approach we take here. In particular in the 
study of a n.n.n. Hubbard model on a square lattice 
as a prototypical model for underdoped cuprates, Fu- 
rukawa et al.Q, and Honerkamp et al.'3|, found a mu- 
tual reinforcement of antiferromagnetic (AF) and d-wave 
pairing tendencies when the Fermi energy was close to 
the van Hove singularities at the saddle points (SP) at 
(±7r, 0), (0, ±7r). They pointed out that the same mech- 
anism is at work in the RG flows in the half-filled 2-leg 
Hubbard ladder (2LHL). At the heart of the mutual rein- 
forcement is the geometrical fact that in these cases a sig- 
nificant fraction of particle pairs that experience strong 
Umklapp scattering in the (tt, tt) channel, e.g. two par- 
ticles close to the SP (tt, O), have small total momentum 
and thus couple into the Cooper channel as well. Since 
both the AF (tt, 7r)-scattering and the large momentum 
transfer pair scattering in the (i2.2_j,2-channel are repul- 
sive, a mutual reinforcement occurs. In the half-filled 
2LHL, this works down to arbitrarily low scales. In 2D, 
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the overlap between the channels decreases for low scales 
such that a finite interaction strength is needed for this 
mechanism. The strong coupling behavior of the ladder 
has been fully analyzed using Density Matrix Renormal- 
ization Group (DMRG)(^ and by bosonization by Lin 
et al.0 who named the groundstate a d-Mott phase. It 
shows the unusual behavior of enhanced correlations si- 
multaneously in the AF, d-wave pairing and d-density 
wave (or orbital AF) channels and is intimately related to 
resonating- valence bond (dRVB) states at large U . Note 
all correlations remain strictly short range due to the 
presence of both spin and charge gaps. Clearly this is a 
case of a nonperturbative groundstate which cannot be 
even qualitatively described by a mean field approxima- 
tion based on long range order. 

We analyze first the case of the 2LHL which is a good 
test case since the groundstate is highly nontrivial and 
turn then to the n.n.n. Hubbard model on a square lat- 
tice. The limitation of Lanczos algorithms to a relatively 
small number of fc-points means that in 2D a representa- 
tion of the whole FS is impossible. However in the RG 
flows of the 2D n.n.n. Hubbard model, the dominant di- 
vergent scattering processes appear in the vicinity of the 
SPs. This allows us to restrict our reduced model, de- 
rived from the RG flow, to processes connecting a discrete 
set of k points located near the two SPs. If this discrete 
set of /c-points is restricted to just the two SPs, then the 
RG flow reduces to a small number of coupled differen- 
tial equations. Two different fixed points were found by 
Dzvaloshinskii [ifij and by Lederer et al.jTJI. The latter 
found a divergent flow to strong coupling when the next 
to leading order particle-hole processes are included. The 
2D RG flows have a similar form for the dominant pro- 
cesses connecting the SP regions. Our numerical results 
for this case reveal both spin and charge gaps at the sad- 
dle points and show enhanced correlations in both AF 
and d-wave pairing channels. They confirm the similar- 
ity to the d-Mott phase of the 2LHL, in line with the 
earlier interpretations. 

Let us briefly describe the numerical approach. An ex- 
tended discussion and results for the 2LHL and other ID 
systems will be published elsewhere [1^. Our scheme is 
based on the numerical analysis of the reduced Hamil- 
tonian Hp^ obtained from the RG flow close to the crit- 
ical scale Ac. 77a, acting on states in a A-shell around 
the Fermi points or surface, is discretized using a finite 
number of fc-points (illustrated in Fig. ^ for two Fermi 
points). Generically it has the form 

k.(T 

where e(fc) denotes the kinetic energy, A is a global cou- 



pling constant, Q. the total volume, VA(fci, fe, ^a) the dis- 
cretized coupling function and fc4 = fci + fc2 ~ ^3 (modulo 
umklapp) due to momentum conservation. The momen- 
tum dependence of the interaction function VA(fci, fc2, fcs) 
is determined by the RG couplings at scale A: 

VA(fci,fc2,fc3) = .9[A] (^Patch(fci),Patch(fc2),Patch(fc3)^ . 

Here g[K\ denotes the ratios of the diverging couplings 
and Patch(fc) is the patch index of fc-point fc. Note 
that in the reduced BCS Hamiltonian, only interactions 
Vj^{ki^k2,kj,) with fc -I- fc2 = are kept. In contrast to 
this, here in order to describe the short-range d-Mott 
state, a whole set of scattering processes with nonzero 
fc -|- fc2 « must be kept in 77a- An analogous width for 
^(^1,^2,^3) for momentum transfers « (tt, tt) is nec- 
essary. To incorporate this essential feature we assume 
that all scattering amplitudes with the same patch in- 
dices take the same value. The global coupling constant 
A is in principle determined by the the initial conditions 
and the final scale A. Here however we use it as a param- 
eter in order to investigate the interplay between kinetic 
energy and scattering processes and to analyze the effects 
of the finite number of fc-points. 

In the next step the discretized Hamiltonian is diag- 
onalized using a Lanczos algorithm. Due to momen- 
tum conservation (modulo umklapp scattering) a num- 
ber Nk ^ 20 of orbitals can be treated. In this way en- 
ergy gaps, static and dynamical correlation functions are 
easily accessible, similar to standard real space diagonal- 
ization algorithms. As a future development the Lanczos 
scheme might be replaced by a DMRG algorithm, thus 
possibly allowing more fc-points for a higher resolution 
on the FS. 

We illustrate the method with an application to the 
two leg Hubbard ladder at half filling. We choose tj_ < 
2f|| in order to have 4 Fermi points. The RG equations 
have been derived and analyzed in j^. In particular for 
general repulsive initial interactions the couplings flow 
towards a fixed ray, i.e. a reduced Hamiltonian with fixed 
ratios among the diverging scattering vertices. 

We characterize the groundstate of the half-filled 2LHL 
first by calculating different gaps: the spin gap: Aspi„ = 
EQ{Ne, S"' = 1) -£;o(A^e, 5"^ = 0), the single particle gap: 
Aip = [Eq{N, - 1, 1/2) + Eo{N, + 1, l/2)]/2 - E^iNe, 0) 
and the two particle gap: Aap = [EQ{Ne — 2,Q) + EQ{Nf.+ 
2, 0)]/2 — EQ{Ne, 0). Results are shown in the left panel 
of Fig. 121 for different system sizes and choices of A. One 
nicely sees that all gaps remain finite, even for — *■ 00. 
Interestingly the spin gap and the two particle gap are ex- 
actly degenerate, while the single particle gap is slightly 
larger. This can be understood as an approximate real- 
ization of the 5*0(8) symmetry of the fixed ray, discov- 
ered in Ref. Sj. In our discretized scheme only a SO (5) 
symmetry seems to be manifest. 
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FIG. 1: Discretization of the A-shell around two Fermi points. 
Small circles denote individual fc-points. A set of fc-points 
belonging to the same branch is called a patch. 
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FIG. 2: Left panel: Spin, two particle and single particle 
gaps in the d-Mott phase of the two-leg Hubbard ladder as a 
function of the interaction parameter A. Right panel: Order 
parameter susceptibilities for various kinds of order in the d- 
Mott phase of the two-leg Hubbard ladder. No evidence for 
long range order is found. However the staggered spin, the 
d-wave pairing and the staggered flux correlations are equally 
enhanced with respect to the noninteracting ground state. 



Now we consider the correlation functions. We inves- 
tigate particle-hole and particle-particle response with 
both s-wave and c?-wave formfactors. The Q — (tt, tt) 
charge density wave structure factors are defined as fol- 



lows: 
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where 



f{k) = 1 for the s-wave channel and f{k) — exp{iky) 
for the d-wave channel. The structure factors for the 
spin density wave and the pairing correlations are de- 
fined similarly. The results for the case X/A — 2 are 
shown in the right panel of Fig. [SJ normalized to the re- 
spective expectation value in the noninteracting Nj: = 8 
groundstate. The response for the standard (tt, tt) spin 
density wave, the (tt, tt) d-charge density wave (orbital 
AF) and the d— wave pairing correlations are equally en- 
hanced. However the finite size scaling does not support 
long range order, in complete agreement with the results 
obtained by bosonization |^. 

The iV-patch RG analysis for the n.n.n. 2D Hubbard 
model revealed a saddle point regime Q where the FS 
is close to the SPs but sufficiently curved such that the 
Brillouin zone diagonals are not nested. Then the leading 
fiow is given by the scattering processes between the two 
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FIG. 3: (a) The four scattering vertices of the two patch 
model, (b) The location of the fe-points for the 8 and 16 
fc-point meshes. The curved lines denote the arcs. 

SP regions. The rest of the FS, which we will call arcs, 
contributes less to the flow and will be considered later. 
First let us focus exclusively on the SP regions. Then 
the flow between the SPs is the one described by Lederer 
et al.^3l and Furukawa et al.0. There are four different 
types of processes, gi, . . . 54, depicted in Fig. O (a). In 
principle the scattering vertex VA(fci, fc2, ^3) depends on 
the precise location of the wave vectors ki in the SP re- 
gions. As a simplification we assume that they have the 
same value. To a certain extent this can be justified by 
the N-patch RG. Note however that in the one- loop flow 
the couplings change smoothly from arc to SP regions. 

We implement the two patch model using the two 
discretizations shown in Fig. |3| (b). From the RG 
approach 7J the ratios of the couplings close to the critical 
scale are known: : 172 : <73 : <74 = : 1 : 2.2 : —1. Using 
this we calculate energy gaps and correlation functions in 
the same way as for the 2LHL case discussed above. The 
gap structure and structure factors are shown in Fig. 0] 
The sSDW and dSC channels are enhanced in agreement 
with the divergences found in the 1-loop RG flow 01 and 
the dCDW is also enhanced. The results show a striking 
analogy to the d-Mott phase of the half-flUed 2LHL. This 
clearly supports the idea that the initial divergence leads 
to the truncation of the FS with spin and charge gaps 
opening near the SP. We examine next the behavior of 
the remaining arcs of the non-nested FS as the RG scale 
is lowered further. These arcs couple to the SP regions 
in the Cooper channel. 

Although there is a charge gap for particle pairs Sp'^s, 
S-p gi in the SP regions, quasiparticle pairs with zero 
total momentum a^r^, a_-j, ^, from the arcs can be scat- 
tered there as virtual excitations. The coupling can in 
principle be derived from the iV-patch RG, its action is 
{k = (iw, k)) 
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FIG. 4: Left panel: Spin, two particle and single particle 
gaps in the two patch model as a function of the interaction 
parameter A. Right panel: Order parameter susceptibilities 
for various kinds of order in the two patch model. No evidence 
for long range order is found. However the staggered spin, the 
d-wave pairing and the staggered flux correlations are equally 
enhanced with respect to the noninteracting ground state. 



Formally integrating out the SP regions gives rise to an 
induced pairing interaetion on the arcs 



^AA 1^ 2^9k(^k,sa^ 



k,s' 9flf^~k' ^s' ^k' ,s • 



Here, Xs the static d-wave pairing susceptibility per 
volume of the SP regions. Xs remains finite at all tem- 
peratures T, as the pair correlations in the d-Mott state 
decay exponentially. However the induced attraction en- 
hances the weak d-wave attraction within the arcs, and 
eventually, at low enough temperatures the arcs will un- 
dergo a superconducting transition. This mechanism re- 
sembles that proposed by Geshkenbein et al. ^3 with 
preformed pairs at the saddle points. As the SP regions 
are insulating, the superfluid weight will be that of the 
arcs alone. With a very similar reasoning for the (vr, tt)- 
processes in the spin channel we can understand how the 
system develops AF order when the arcs become more 
nested at sufficiently small doping. 

We may further ask what criterion determines whether 
the superfluid density is reduced in the groundstate. The 
answer lies in the strength of the mutual reinforcement 
of Cooper and (tt, 7r)-channels at the SPs. If the (tt, tt)- 
'almost-Cooper' processes enhance the d-wave pair scat- 
tering at the SPs more than the coupling in the Cooper 
channel to the arcs, we can expect that the system fa- 
vors the reduced superfluid density. In this situation the 
energy gain through the mutual reinforcement - manifest 
in spin and charge gap - is larger than the loss of su- 
perconducting condensation energy by the partial gap to 
Cooper pairs in the SP regions. If the mutual reinforce- 
ment is not strong enough compared to other processes. 



a d-Mott or equivalently dRVB condensate will not ap- 
pear as the initial instability. Instead a long range order 
such as ordinary d-wave superconductivity on the whole 
FS at larger doping, when the FS moves away from flat 
SP regions, or AF order at small doping - sacrificing the 
partial spin gap - when the nesting over the the whole 
FS becomes too strong. 

Summarizing we have determined a reduced fermionic 
Hamiltonian that contains the essential ingredients of a d- 
Mott state. For the half-filled two-leg Hubbard ladder we 
have shown that exact diagonalization of this Hamilto- 
nian reproduces the known spin and charge gaps. For the 
2D n.n.n. Hubbard model we have described how such 
a Hamiltonian captures the leading scattering processes 
between the saddle points, and hence how a d-Mott con- 
densate can truncate the Fermi surface near the saddle 
points. Finally we have outlined how superconducting 
and antiferromagnetic states can arise as additional in- 
stabilities on the remaining parts of the Fermi surface, so 
that the d-RVB phase can be viewed as an unstable fixed 
point as proposed by Anderson 0|. In 2D our picture 
still relies on several approximations and assumptions, 
but we hope that the present work may provide some 
guidance for further research. 
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